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Abstract 

We present a new calculation of the 0{a 2 s ) coefficient of the energy-energy correlation 
function (EEC) using two different schemes to cancel infrared and collinear poles. The 
numerical evaluation uses the phase space slicing and the hybrid subtraction method. 
Both schemes converge with decreasing slicing cut. The results are independent of the 
scheme for small cuts. For the pure phase space slicing method, the cut must be below 
10 -6 to achieve good results. All four approaches agree with each other and confirm the 
results of Kunszt and Nason and Glover and Sutton, for the latter also with respect to 
contributions of different colour factors. 
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1 Introduction 



The energy-energy correlation function E has been used by all for LEP experiments [pj and the 
SLD experiment at SLAC to measure the strong coupling constant a s in e + e~ annihilation 
at the Z resonance. £ is defined as a function of the angle x between two particles % and j in 
the following form 

dcOSX AcOSX^ent jvtltife ^ 

where and Ej are the energies of the particles and E is the total energy of the event, E 2 = s. 
The sum runs over all pairs i, j with cos \ m a bin of width A cos x- cos X ~ A cos x/2 < cos x < 
cosx + Acosx/2. Each pair enters twice in the sum. The limits Acosx — > and N event —> oo 
have to be taken in (1). a is the total cross section for e + e~ — > hadrons. 

In perturbative QCD the energy-energy correlation function (EEC) is given as a series in 
a s which we write as 

1 dTj a s (/i) ( a s (jj)\ 2 f ( ji 



aoimX 2, A ^ + {^) ^{%)«x) + B M ) + Otf), (2) 

where fa = {11N C - 4T R )/3. For QCD we have C F = 4/3, N c = 3 and T R = N f /2 where N f 
is the number of active flavours at energy E. 

The first order term A(x) is calculated from the well-known one-gluon emission diagrams 
7*, Z — > qqg, where x is th e angle between any of the three partons. It has been first calculated 
by Basham et al. with the result 

A( x ) = C F (1 + ^) 3 ^^ ((2 - 6a; 2 ) ln(l + l/u) + - 3) (3) 

where u = cot 2 x/2. The next-to-leading order (NLO) contribution B(x) is obtained from the 
processes 7*, Z — > qqg at one loop and 7*, Z — > qqgg, qqqq at tree level. Several groups have 
calculated B(x) i n the past: Richards, Stirling and Ellis (RSE) 0], Ali and Barreiro (AB) ||, 
Schneider, Kramer and Schierholz (SKS) 0, Falck and Kramer (FK) |7]], Kunszt and Nason 
(KN) || , Glover and Sutton (GS) § and just recently Clay and Ellis (CE) |K|. Most of the 



calculations, except those of SKS and FK, are based on the ERT matrix elements |irTJ. The 
numerical evaluations differ by the method used to cancel infrared and collinear divergences. 

Considering the results for B(x) obtained by RSE, AB, FK and KN one notices that they 
differ at the special point x — 7r /2 by roughly 50%. The SLD collaboration used all four 
theoretical calculations to determine a s from their measurements. They averaged the values of 
a s obtained this way and increased the theoretical error accordingly. This additional error of a s 
due to the uncertainty over which the coefficient B(x) in the EEC is correct is equal to ±0.006, 
which is appreciable, considering that the pure experimental error is only (+0.002, —0.003) and 
the error from the hadronisation corrections is ±0.002. This unsatisfactory situation has raised 
renewed interest in the calculation of the EEC function and Glover and Sutton || have pre- 
sented a new calculation of B(x) using three methods, called subtraction, phase space slicing 
and hybrid subtraction method. They confirmed the result of Kunszt and Nason H. In the 



most recent calculation of Clay and Ellis |T0[ a larger value of B(x) is obtained which agrees 
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with the result of Falck and Kramer 0, but is in disagreement with the other recent calcu- 
lations of GS and KN. CE argued, on the basis of their much higher statistical precision of 
±0.3%, that the theoretical error on a s could be reduced now; however, this does not seem 
justified since the spread of theoretical predictions did not change with this new calculation. 
One might think, that the theoretical prediction of the 0(a 2 s ) coefficient B(x) should be unique 
with no other uncertainty than the error of the numerical evaluation. However, we we will see 
later that this need not always be the the case. 



In order to compute the NLO term B(x) in (2) it is necessary to combine the contribu- 
tion from the 3-parton one-loop diagrams (multiplied with the LO graphs) with the 4-parton 
7*,Z — > qqgg,qqqq processes. The virtual matrix elements contain infrared and collinear sin- 
gularities which cancel with the singularities in the matrix elements of the 4-parton final state. 
This cancellation of the infrared and collinear poles is done analytically using dimensional reg- 
ularisation. The extraction of the infrared and collinear poles in the regulator e (dimension 
n = 4 — 2e) and the complete calculation of the finite terms for e — > has been done by 
two groups independently with identical results [IT, [IS]. An analytical calculation of the full 



4-parton expressions at finite e is not possible. Therefore these expressions are simplified in 
such a way that they contain all the infrared and collinear singularities. Analytical integrals 
of these simplified 4-parton matrix elements were then added to the virtual 3-parton contri- 
butions and the sum was shown to be finite in the limit e — > 0. The difference between the 
exact 4-parton matrix elements and the simplified expressions was then computed in 4 dimen- 
sions using numerical methods. To calculate this difference, essentially two methods have been 
applied which are denoted subtraction and phase space slicing method in the literature fTi 



These two methods can nicely be explained following Kunszt and Soper [TJ]] by considering a 
simple one-dimensional integral 

T — lim ( C —x e F(x) - -F(0) ) (4) 
<^o \Jo x e J 

where F(x) is a known function representing the 4-parton matrix elements, x is the variable 
which produces the singularity, i.e. it is the energy of the gluon, the angle of two partons or 
an invariant mass which approaches zero in the infrared or collinear limit. The integration 
over x represents the additional phase space of the parton which produces the pole term. The 
integrand is regularized by the factor x € as it appears in dimensional regularization. The first 
term is still divergent for e — ► 0. So to cancel the divergence of the second term, which represents 
the divergent contribution from loop diagrams for the 3-parton matrix elements, one must do 
the integral for e ^ 0. Now the two methods to perform this integration are as follows. First in 
both methods one isolates the singularity of the integrand by subtracting the pure pole term 
F(0)/x and adding it to the second term with the result 

T = lim [ ( l —x e (F(x) - F(0)) + F(0) C —x e - -F(0)} 

e -^o \Jo x Jo x e J 

= f X -(F(x) - F(0)) (5) 
Jo x 



This way the infrared singularity is cancelled and one is left with a manifestly finite integral 
which is evaluated at e = 0. Due to the complicated structure of the 4-parton matrix elements 
this last integral must be calculated numerically. The above procedure defines the subtraction 
method and was first used for the calculation of event shape distributions in e + e~ annihilation 



3 



by Ellis, Ross and Terrano 
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An alternative approach is the phase space slicing method which is familiar from many 
QED calculations. It was first applied in a NLO calculation of the thrust distribution in e + e~ 




annihilation by Fabricius, Schmidt, Schierholz and Kramer [I5|]. In this method the integration 
region in (4) is divided into two parts, < x < y min and y min < x < 1. In the first region, the 
function F(x) can be approximated by F(0) provided the arbitrary slicing parameter y m i n <C 1, 
so that 

T — lim I / — x e F(x) + / — x e {F(x) - F(0)) + F(0) / — x e - 

\Jymin x Jo x Jo X 

= / —F(x) + / — - F(0)) + F(0) lny min 

^ Drain ** ^ 

~ ^ F ( x ) + F (o)l n2/min . (6) 

^ ymin 

The last equation in (67) follows when y m in < 1 is chosen small enough, so that the finite 
integral 

rVmin (],X 

T f = — (F(x) - F(0)) (7) 

Jo x 

can be neglected. If this is the case the integral T should not depend on y min which usually 
is taken as evidence that Tf is indeed negligible. This cutoff dependence was investigated, for 
example, for the thrust distribution in |T6[ and for y min small enough (y m in < 10 -4 , where 



ymin was the invariant mass cutoff) the result was found to agree with the thrust distribution 
obtained with the subtraction method ]T7| I8|. It is clear that if the finite contribution Tf in 



(7) is not neglected both the subtraction and the phase space slicing method must give the 
same results, independent of how small the cutoff y m in is. If the finite term is kept, we have for 
T 

/l Ax 
—F(x) + F(0)\ny min + T f . (8) 
.y'min 

This version is called the hybrid subtraction method in the recent work of Glover and Sutton 
P| and will be called the hybrid method in the following. 



It is clear that in this simple example, the results for T in the subtraction method and in 
the hybrid method are identical, independent of how y min has been chosen in (8). However, 
in the actual application to the EEC the situation is more complicated. In this case the 
function "F(x)/x" is the matrix element of the process e + e~ — > 4 partons which depends 
on five variables, whereas "F(0)" corresponds to a contribution to the 3-parton final state, 
e + e~ — > 3 partons, which depends only on two variables. Both contributions have to be 
integrated out up to the one variable cos \ and in order to complete the definition of the EEC 
function one has to supply a prescription for the calculation of cos %■ m t ne region x > y min 
where x stands for one of the variables for which the matrix elements become singular, it is 
most natural to determine cosx from the pairing of two of the 4 parton momenta Pi,P2,P3 and 
p 4 , i.e. 

cos x = PiPj (9) 
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with i,j = 1,...,4. The finite integral Tf in (7) can be calculated at least in two ways: (i) 
the integrand "F(x)/x" is calculated with 4-parton kinematics, i.e. with (9) and i,j = 1, ...,4, 
whereas the subtracted integrand "F(0)/x" is calculated with variables corresponding to 3 
partons in the final state, i.e. with (9) and i, j = 1,..,3, of course after the integration over 
x and two more variables is performed in the 4-parton case; (ii) the whole integral with inte- 
grand (F(x) — F(0))/x is evaluated with 3-parton kinematics. The latter procedure amounts 
to performing the integration over x and two more of the 4-parton variables with two variables 
considered constant which are then identified with the two variables describing a 3-parton final 
state. With other words, for the evaluation of Tf one has the freedom to treat the contribution 
F(x)/x exactly or with recombination of two of the partons into one jet. The selection of the 
partons i,j, which are recombined is determined by that one of the variables yij = (pi +Pj) 2 /s 
which leads to the singularity and which is identified with x in the example integral. It is clear 
that this recombination procedure is not unique and we shall consider two possibilities later on. 
In the limit y m in — > the distinction (i) and (ii) should not matter. But for finite y m in > we 
expect a difference between (i) and (ii). The procedure (i) is used in the subtraction approach, 
so that for (i) the result of the hybrid method should be independent of y m in whereas with (ii) 
the region where the recombination is performed is changed so that we expect that the total 
integral will depend on y m i n . The y min dependence should be smaller than in the phase space 
slicing approach, in which the integral Tf is neglected altogether. It is the general consensus 
that all these methods should give legitimate results, provided y m in is chosen small enough. 



As in previous work we focus on the zero- resolution limit (y m in — > 0), i.e. we do not apply a 
recombination procedure to all of the two-parton combinations to form jets before the calcula- 
tion of the EEC function is done. Calculations with a jet recombination for all possible parton 
pairings were performed by AB [|] and by SKS [[J using Sterman- Weinberg (e, 5) cuts and also 
in a recent experimental analysis of the ALEPH collaboration 119 



To shed some further light on the question where the discrepancy between the various 
calculations could arise we have made the effort to perform a new calculation of B(x) using 
the hybrid method in the form (ii) described above. In recent work |2(| we have studied the 



asymptotic behaviour of B(x) for % — > 7r in order to compare with the predictions of the leading 



logarithm approximation for large angles [21] using the hybrid method in the same form. Here 
we extend these calculations to the whole x range. In addition, we separated the contributions 
of the different phase space regions, which can be characterized by the first, second and third 
term in (8) for the example integral. Furthermore to explore the source of the discrepancies we 
decompose B(x) into the contributions from different colour factors 

B( X ) = C 2 F BcAx) + C F N c B Nc ( X ) + C F N f B Nf ( X ) (10) 

which can be compared with the results of other calculations. In section 2 we describe the two 
methods used to calculate the 0(a 2 s ) contribution to the EEC and present the results for the 
various colour factors separated according to the different phase space regions. In section 3 we 
compare our results with the earlier calculations and comment on the different approaches. We 
end with some concluding remarks. 
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2 Calculational Methods and Results 



The new calculation of the EEC in O(a^) is based on the approach described in detail in It 



starts from known 0(a 2 s ) matrix elements for e + e — > 3partons [|TT], [13| and e + e — > Apartons 



23[ pTTfl . To obtain finite results, in which all infrared and collinear singularities are cancelled 
we introduce a phase space slicing cut to separate the 4-parton phase space from the region 
in which the integration over one of the invariants producing the singularity has been done 
analytically in n dimensions. This integration region is defined by 

Vij = (Pi + Pjf/S < Umin (11) 

where y m in is the parameter to separate the two regions and i, j are the labels for two of the 
4 parton momenta pi,i = 1,2,3,4. This slicing procedure with an invariant mass cut is most 
convenient for analytical calculations. The integration over only one of the invariants is possible 
only when the singular contributions have been separated. To achieve this a partial fraction 



decomposition of the 4-parton matrix elements is carried out (see [22| for details). Then the 
EEC function is calculated from three separate contributions: 

(a) The first contribution contains the singular parts of the partial fractioned 4-parton matrix 
elements integrated over one invariant inside the slicing cut y m in together with the virtual 
corrections to the qqg final state. This contribution corresponds to the term F(0) \ny m i n in (8) 
for the example integral and depends strongly on y m in with the dominant term oc (—In 2 y m i n ). 

(b) The second contribution contains the remaining non-singular part of the 4-parton matrix 
elements which is the difference between the full 4-parton matrix element and the singular 
part already included in (a) integrated over the same region as in (a). This part corresponds 
obviously to the term Tf in (8). 

(c) The third part consists of all contributions outside the singular region (a), i.e. with the 
intgration over y^ > y m i n , which is computed numerically with cosx given by (9) with i,j = 
1,2,3,4. This part represents the first term in the example integral on the right-hand side of 
(8). 

The calculation of the 4-parton contributions needed in (a) proceeds as follows. First the cross 
section for e + e~ — > qqgg (the qqqq final state is less singular and does not require partial 
fractioning; otherwise it is treated similarly) has the general form 



a, N 2 



d'o = f( ytj )dPS®. (12) 

The right-hand side of (12) has pole terms proportional to y^ , (ij = 13, 14,23,24,34) which 
are separated by the partial fractioning. For example, the contribution proportional to the 
colour factor C\ has the structure 

fiVij) = — + (1 <-> 2) + (3 <-> 4) + (1 <- 2,3 <- 4) (13) 

2/13 

where the momentum labels are 1,2,3,4 = q,q,g,g. The terms proportional to y^ 1 and y^, 
respectively y^ and y^i, become singular when one of the gluons is infrared or collinear with 
the quark, respectively antiquark. They produce the dominant negative contributions to B(x) 
after integration over the unresolved regions yi 3 < y m i n , 2/23 < Vmin, etc. when they are added 
to the virtual contributions. This integration is done only for one of the four terms in (13) 
which are related by permutation of the momentum labels. It is important to note that with 
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this procedure of slicing we do not calculate genuine 3-jet cross sections since only one of the 
pairings ij of two partons is recombined into one jet in the first, second etc. term in (13). The 
four terms in (13) are separated into singular and non-singular terms. The singular terms are 
regularized by dimensional regularization. The singularities in e after integration compensate 
against the singularities in the 0(a 2 s ) one-loop corrections to e + e~ — > qqg. The result from the 
contribution (a) will be denoted the singular contribution, abbreviated by sing in the figures 
and tables which contain our results. 



The other two contributions (b) and (c) come exclusively from the 4-parton tree diagrams. 
The part (b) consists of the non-singular terms in A 13 /y 13 etc. in (13) which are also integrated 
up to the cutoff. They do not participate in the cancellation of the infrared/collinear singular- 
ities between tree and virtual diagrams. The third part (c) is connected with the contribution 
above the slicing cut, 2/13 > y min (in the first term of (13)), which is computed numerically. 
It is important to note that the same expression for f(yij) with partial fractioning is used as 
in the singular region and that the integration up to y m in over t/i 3 is applied only to the first 
term in (13). In the other terms the integration is over the variables 2/23, 2/i4 and 2/24- The par- 
tial fractioning is not unique in the sense that non-singular terms can be distributed between 
the terms proportional to Ay in (13). In the figures and tables and in the discussion below 
the contributions (b) and (c) will be called finite terms (fin) and real terms (real), respectively. 



To proceed with the calculation we must specify the singular contributions given by A^jy^ 
in (13). These singular terms factorize into the pole terms times a factor which can be 
identified with the LO matrix element for e + e~ — > qqg. Only because of this factorization 
one is able to perform the cancellation with the virtual contributions. The identification of 
this factor with the LO matrix element amounts to specifying combinations of the 4-parton 
kinematic variables for the qqgg final state which define a 3-jet final state. These 3 jets can 
then be identified with the 3 partons of the LO matrix element. Usually this is done with the 
help of invariants built from the momenta of the qqgg final state. For example in the ERT 
approach |lT| which was followed in all the calculations using the subtraction method, the rela- 
tion was such, that in the pole term proportional to y~^ the 3-jet variables are 7/134, 2/24 and 1/123 
(yijk — (Pi +Pj +Pk) 2 / s). For three massless jets one has 7/134 + 2/24 + 2/123 = 1- This agrees with 
the 4-parton energy-momentum relation 7/134 + j/ 2 4 + 2/123 — 2/13 = 1 only in the limit 7/13 — > 0. 
Therefore in [22| two other schemes have been considered: (1) the so-called KL scheme where 
2/134)2/24 and 2/123 — 2/13 are used instead; (2) the so-called KL' scheme, where the 3-jet variables 
are 2/134,2/24 — 2/13 and 2/123- Of course there are many more possibilities. This non-uniqueness 
of the definition of 3-jet variables for a 4-parton final state can not be avoided if one wants 
to cancel the infrared and collinear singularities between real and virtual contributions before 
the distribution in the final observable, i.e., in our case, in cosx, is computed. This difference 
in the 3-jet variables is supposed to have no effect for y min — > 0, i.e. in this limit we expect 
no difference in the results for the KL, KL' and ERT scheme. Since y m in is really never put 
equal to zero we must anticipate an effect from these different schemes in practice. From this 
discussion it is clear that in the calculation of the EEC function there is an ambiguity which 
is not visible in the simplified example integral discussed in the introduction. There, F(0) can 
be defined uniquely: since F depends only on one variable, F(0) is the residue of the pole. 
With more than one variable, however, the equivalent of F(0) depends on which of the other 
variables are kept fixed and identified as 3-jet variables. 
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In the following we present separate results for the contributions of the C F , CpNc and CpNf 
parts of B(x) as defined in (10), where we include the colour factor, i.e. we plot sin 2 x C f B Cf (x) 
etc. for the three contributions sing, fin and real for a slicing cut y m in — 10~ 4 . We have 
chosen the KL' scheme for these plots and factored out sin 2 x to be able to plot with a linear 
scale. In Fig. 1 we have plotted the contributions for the Cp part as a function of cos x- We 
see that the real part is large and positive whereas the singular part is large and negative as we 
expect. The sum of both terms is small and negative for cosx —> ±1, since the singular terms 
oc ln n ymin, n — 1; 2, cancel in the sum. The finite term (fin) is very small. In the figure it is 
hard to distinguish it from the zero line. The same plot for the CpNc part is seen in Fig. 2. 
Here we have two classes of pole terms coming either from the qg{qg) or the gg recombination. 
The corresponding real and non-singular parts are denoted real (13), f in(13) and real (34), 
fin (34), respectively. In the singular term both contributions (13) and (34) are added together 
for convenience. We see that the fin (34) is larger than fin (13) and real (34) (real (13)) is 
large positive (large negative). The total sum is positive. For the CfNj contribution, calcu- 
lated for Nf = 5 and exhibited in Fig. 3, we have only the singular and the real contributions. 
The finite terms are included already in the singular terms (see for details). The singular 
and the real contributions contain terms which are proportional to (—lny min ) and \ny min , re- 
spectively. These terms cancel in the sum which is negative and smaller. In addition we have 
plotted the contribution of the qqqq interference term. This term is actually proportional to 
Cf(Cf — Nc/2) and should be distributed to the Cp and CfNc terms. For convenience we 
have included it here. Since it has no singularities it can be calculated numerically without 
difficulties. Its contribution is very small and negligible compared to the other terms. The sum 
in Fig. 3 includes this interference term, i.e. the sum contains all 4-quark contributions. In Fig. 
4 the C F) CpN c and the CpNf contributions are plotted together with the sum of all three 
colour terms. We see that the CfNc term is the most important one and it is positive. The 
CfNj term is smaller and negative. The Cp term is less important. The sum is essentially 
given by the sum of the CpNc and the CpNf parts except near the wings of the EEC distri- 
bution. The corresponding curves for the KL scheme look similar and will not be shown here. 
For larger y m i n values the results change. For the KL' scheme the EEC correlation function 
decreases with increasing y m in- For the KL scheme this decrease occurs in particular for the 
wings of the NLO EEC function. 

To obtain an overview of the most interesting central region we give the numbers in Tab. 
1 for cosx = and y min = 10~ k , k = 2,3,4,5,6. Since the numerical results fluctuate from 
bin to bin we have made a polynomial fit to our results in the interval —0.40 < cosx < 0.40 
and quote the result of this fit at cosx = 0. We give the results for the different colour factors. 
We see that the CpNf term is independent of y m in, which we expect because of the less sin- 
gular behaviour of this part. The variation of the EEC function results in particular from the 
variation of the C\ term, whereas the variation of the CpNc contribution with y m i n is much 
smaller. In the KL case the variation of the Cp contribution is smaller, so that also the total 
sum is fairly constant. The splitting up of B(x) into the contributions from the various regions 
for cosx = and y m in = 10~ 4 is given in Tab. 2 for the KL' case. We see that the finite part 
fin (34) is still not negligible for this slicing cut. It is still more than 10% of the total sum and 
larger than the absolute value of the total Cp contribution. The dependence of this particular 
contribution on the value of the slicing cut is shown in Fig. 5. The qqqq interference term is 
negligible and much smaller than the total error. 
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Results for the pure phase space slicing method can be obtained by subtracting the contri- 
butions of the finite terms from the results in Tab. 1. The resulting numbers for the same y m i n 
as in Tab. 1 are collected in Tab. 3, again for the different colour factors and the KL' and KL 
schemes. We see that the convergence with y m in is slower than for the hybrid method. The 
slow convergence comes mostly from neglecting the term f in(34) which is the largest one of 
all the finite contributions (see Tab. 2). Even for the larger y m i n values the results for the KL' 
and KL approaches differ very little. The difference between the two schemes in Tab. 1 comes 
essentially from the finite term in the C\ contribution which changes stronger with y m i n in the 
KL' scheme. The slower convergence of the pure phase space slicing method was also observed 
by Glover and Sutton ||. If we compare the results for the hybrid and the phase space slicing 
method at y m in — 10~ 6 we see that they agree inside the numerical errors. The difference 
of B(x — 7r /2), however, is still 1.15. This comes essentially from the term fin (34) which 
decreases very slowly with decreasing y m i n , but is still non-negligible. In the calculation we 
obtained f in(34) = 1.154 ±0.006 for both the KL' and the KL case. So, even with a slicing cut 
as small as 10~ 6 the finite parts are still not very small and produce an error of approximately 
2% for B(x — 7r /2) if one uses the pure phase space slicing method. At y m in = 10~ 5 the B(x) 
differ already outside the given error. Thus, to obtain good values for the NLO coefficient B(x) 
with the phase space slicing method one must go to extremely small y m in of the order of 10~ 6 
or smaller. The dependence of B(x = tt/2) on y m in is shown in Fig. 6 for the hybrid and space 
space slicing method and for the KL' and KL scheme separately for y min between 10~ 6 and 
10 -2 . We observe that there is little difference between the KL' and the KL scheme in the case 
of the phase space slicing procedure since these two schemes differ essentially only in the finite 
parts. For the phase space slicing method, B(x = tt/2) increases monotonically with decreasing 
ymin- With the hybrid method B(x) first increases to a maximum near y m i n — 10~ 3 and then 
decreases towards the final value at y m in = 10~ 6 . For y min < 10~ 4 the two schemes, KL' and 
KL, give the same results inside the errors. We consider the value at y m in = 10 -6 which is 
B(x = tt/2) = 50.1 ± 0.9 as our final value. In Fig. 7 we have plotted B(x) as a function of x 
in the range —0.96 < cos x < 0.96 for the same y min for KL', as well as the result for y min = 10 -4 . 



3 Comparison with other Results 

In this section we compare with the results of previous calculations. Such a comparison was 
done already by Glover and Sutton who compared their results with those of AB, RSE, FK 
and KN. Therefore we can restrict ourselves to a comparison with the most recent evaluations, 
namely those of KN, GS and CE. It turns out that our results agree with those of KN and GS 
but not with the results of CE. Kunszt and Nason || calculated B(x) with the subtraction 
method by reorganizing the ERT matrix elements to give numerically stable results. Glover 
and Sutton || determined B(x) with all three methods, subtraction, phase space slicing and 
hybrid method. For the evaluation with the subtraction method they used the ERT matrix 
elements as given in the original work jy]. For the phase space slicing and the hybrid method 



they have constructed a completely independent program based on 2~4"j but using squared ma- 



trix elements rather than helicity amplitudes as in p4| . Our calculation is based on the work 
reported in [^IJ which was developed independently of all the earlier results and which relied 
on the introduction of a slicing cut to isolate the infrared and collinear divergences. A further 
important ingredient was the use of partial fractioning in all invariants yij which lead to sin- 
gularities so that only for one of these variables a slicing cut had to be introduced. We believe 



9 



that a similar technique was employed by Kunszt and Nason and by Glover and Sutton in their 
calculations. 



In Fig. 6 we can see that at larger y min the results vary with i/ mi „, independently of the 
method that was used; even for y m in < 10~ 5 we still have a non-negligible variation. Since 
on the other hand Glover and Sutton give explicit numerical results only for y m i n = 10~ 5 we 
have selected our results for the same y m in to perform the comparison. Since our calculational 
methods and those of Glover and Sutton are not the same, there is actually no reason that the 
two results should be identical for this particular y m i n ; however, we believe that a comparison 
at the same y m in is more sensible than performing the comparison at different values. Since 
Kunszt and Nason used the subtraction method, there is no cut dependence in their results. 

The comparison of our results with the results of these two groups is collected in Table 4. 
We see that we have perfect agreement. In particular our hybrid results agree very well with 
those of Glover and Sutton inside their respective errors (compare KS^ 5 \KS^ with GS^). 
Concerning the phase space slicing method there is a small difference with GS outside the errors 
(compare KS^ 7 \KS^ with GS^). Knowing that the finite parts are still non- negligible at 
Vmin = 10~ 5 , we can not expect perfect agreement in this case. 

Another important check between the various results is the decomposition of B(x) into 
contributions for different colour factors. In Tables 1 and 3 we have given these contributions 
already for B(x) at x = 7r/2. Equivalent results are not available from the publications in refs. 
|8], ||. After their publication, Glover and Sutton calculated the contributions for the different 
colour factors for a comparison with the results of Clay and Ellis. Their results for B(x = tt/2) 
25| are: 



C 2 F B Cf = -2.42 ± 0.92, (14) 
C F N c B Nc = 77.20 ± 2.08, 
C F N f B Nf = -24.97 ± 0.53. 

The sum of these terms is 49.8 ± 2.3. By comparing with our results in Table 1 or Table 3 
we observe perfect agreement of our colour decomposition with the Glover and Sutton results 
given above. 



Concerning the more recent results of Clay and Ellis |H| there is no agreement with our 
results and with those of Kunszt-Nason or Glover-Sutton, respectively. From Fig. 1 of the 
Clay-Ellis paper we read off a value of about 57.0 for B(x = tt/2). The actual precision of their 
numerical results is 0.3% as stated by the authors. So, our value of B as well as that of KN or 
of GS differ by approximately 15% which can not be explained by numerical uncertainties in 
any one of the calculations. In our case the error is below 2% (for y m in = 10~ 6 ). The difference 
is in the contributions with the colour factors Cp and CpN c , since concerning B N Clay and 
Ellis and Glover and Sutton achieved agreement jnj. It is clear that the CpNf term is less 



problematic since the contributions to this term are less singular than those to the Cp and 
CfNc terms. The origin of the disagreement is unclear. Unfortunately we could not study this 
further, since Clay and Ellis have not explained the details of their calculation. 

With respect to the earlier calculations of AB || , RSE || and FK |7J we made some effort 
to understand the differences with these earlier results. The calculations of FK were done with 
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the phase space slicing method with y m i n = 10~ 4 and with the important difference that the 
real contributions are calculated with cuts y^ > applied to all invariants of the 4-parton 
momenta. This was necessary since for the real contributions partial fractioning as in (13) was 
not performed. We found out, however, that for y min = 10~ 4 these additional cuts have a non- 
negligible effect when we apply them to the partial fractioned real term as given, for example, 
in (13). In the present calculation we used cuts only for those variables which are singular in 
the first, second etc. terms in (13). In |7| the reduction of the non-singular terms through 
these additional cuts was compensated by additional terms in the singular contributions, but 
presumably not fully. In the moment it is unclear, whether with much smaller cut values than 
10 -4 the method used by FK would give the same results as the present calculation. From the 
results given by FK it is clear that y min = 10~ 4 is not sufficiently small to obtain correct results. 
Furthermore, from Tables 1 and 3 we see that at y m in = 10™ 4 there is already a difference of 
AB(x — 7r /2) = 6.03 and 6.45 for the KL' and KL schemes, respectively, between the hybrid 
and the phase space slicing method. The results of AB and RSE Q are obtained with the 
subtraction method and therefore should agree with the results of KN and GS using the same 
method. We suspect that in these older calculations the so-called finite pieces (see (5)) are 
calculated also with cuts on all invariants yij in the 4-parton term (first term in (5)). In the 
work of Ali and Barreiro (first reference of [[J) such a cut y^ > 10~ 7 is explicitly stated and 
this cut was also used in the second paper |5], [26|]. Concerning RSE we have no information 



on this point. If the interpretation that the additional cuts in AB and RSE are responsible for 
the reduction of their results is correct, it would mean that, similar as we observed for the FK 
calculation, one would need even smaller cuts below 1CT 7 to achieve convergence. 



We also compared the EEC coefficient B(x) for the other angles as shown in Fig. 7 with 
the results of Kunszt and Nason and found good agreement. We therefore conclude that the 
KN calculation gives the correct NLO coefficient B(x). The older results RSE, AB and FK 
should be disregarded and not be used for a determination of a s . If for the time being we also 



leave aside the results of Clay and Ellis [[K]] since their result has not been confirmed by other 
calculations, the theoretical error on a s mentioned in the introduction can be eliminated. Then 
the SLD analysis || yields the following value for a s from the measurement of the EEC: 

a s {M 2 z ) = 0.125±g:gg| (exp) ± 0.012 (theory). (15) 

This agrees completely with the conclusion of Glover and Sutton M. Now the theoretical error 
comes exclusively from the scale change which in the SLD work originates from a rather 
large variation with a scale factor / in the interval 0.002 < / < 4. 



In conclusion, motivated by the recent results of Clay and Ellis which disagreed with the 
earlier results of KN and GS, we have reevaluated the EEC using a completely independent 
approach which allowed us to obtain results for the phase space slicing and the hybrid method. 
The pure phase space slicing method needs extremely small slicing cuts, smaller than 10~ 6 , in 
order to achieve a good accuracy below 2 %. This is due to the finite terms which converge 
very slowly to zero with decreasing cut. This shows that configurations with extremely small 
invariant masses can contribute to the EEC at all angles, i.e. also at angles far away from the 
two-jet region at x — ft- m our opinion the results of the older calculations seem to suffer from 
a too large slicing cut (FK) or from additional cuts in the evaluation using the subtraction 
method (RSE, AB) and therefore could yield only estimates of B(x). 
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Table Caption 

Tab. 1: NLO coefficient B(x) at cosx = for various y m in and colour contributions in the 
KL' and KL scheme using the hybrid method. 

Tab. 2: Contribution to the NLO coefficient B(x) at cosx = for y m i n — 10 -4 for different 
colour factors from the regions sing, real and fin with the hybrid method and KL' 
scheme. 

Tab. 3: Same as in Table 1 with the phase space slicing method. 

Tab. 4: NLO coefficient B(x) at cosx = for different calculations; (1) Table 3 of 0, (2) 
subtraction method, (3) phase space slicing method, y min = 10~ 5 , (4) hybrid method, 
y m in = 10~ 5 of ref. ||, (5) and (6) phase space slicing method for KL' and KL 
scheme from Table 3, (7) and (8) hybrid method for KL' and KL scheme from Table 

1 {Vrmn = 10~ 5 ). 



Figure Caption 

Fig. 1: Cp contributions to the NLO coefficient sin 2 \ B(x) as a function of cos x for y min = 
10~ 4 in the KL' scheme. Long-dashed line: singular contribution, short-dashed line: 
real contribution, dotted line: finite contribution. Full line: sum of all Cp contribu- 
tions. 

Fig. 2: CpNc contributions to the NLO coefficient sin 2 x B(x) as a function of cosx for 
ymin — 10 -4 in the KL' scheme. Long-dashed line: singular contribution; real con- 
tributions: short-dashed: real (13), dotted: real (34); finite contributions: long 
dashed-dotted: fin(13), short dashed-dotted: fin(34); full line: sum of all CpNc 
contributions. 

Fig. 3: 4-quark contributions to the NLO coefficient sin 2 x B(x) as a function of cosx for 
ymin = 10~ 4 in the KL' scheme. Long-dashed line: singular contribution, dotted line: 
real contribution; short-dashed line: 4-quark interference contribution. Full line: sum 
of all 4-quark contributions. 

Fig. 4: Colour decomposition of the NLO coefficient sin 2 x B(x) as a function of cosx for 
Vmin — 10 -4 in the KL' scheme. Long-dashed line: Cp contribution, short-dashed 
line: CpN c contribution, dotted line: CpNf contribution; full line: sum of all con- 
tributions. 
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Fig. 5: y m i n dependence of C F N c -i in (34) contribution to the NLO coefficient B(x = 7r/2). 

Fig. 6: Umin dependence of the NLO coefficient B(x = vr/2). Upper curves are for the hybrid 
method, lower curves for the pure phase space slicing method. The two curves in 
each set are for the KL' and the KL scheme. 

Fig. 7: The complete NLO coefficient B(x) as a function of cosx obtained with the hybrid 
method in the KL' scheme for y min = 10~ 4 (histogram) and y m in = (points with 
error bars). 
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Table 1: 



NLO Coefficient B( X ) for cos^ = 


Umin 


method 


r 2 


C F N C 


C F N f 


sum 


l(T d 


KL' 
KL 


-2.68 ±0.68 
-2.57 ±0.68 


78.16 ±0.60 
78.03 ± 0.60 


-25.388 ± 0.071 
-25.388 ± 0.071 


50.11 ±0.91 
50.09 ±0.91 


10" 5 


KL' 
KL 


-2.09 ± 0.44 
-2.02 ± 0.44 


79.22 ± 0.36 
79.12 ±0.36 


-25.159 ±0.055 
-25.159 ±0.055 


51.96 ±0.58 
51.87 ±0.58 


10- 4 


KL' 
KL 


-2.18 ±0.26 
-1.83 ±0.26 


80.81 ±0.19 
80.87 ±0.19 


-25.116 ±0.041 
-25.116 ±0.041 


53.52 ±0.32 
53.93 ± 0.32 




KL' 
KL 


-3.77 ±0.13 
-1.61 ±0.13 


81.621 ±0.092 
81.575 ±0.093 


-24.965 ± 0.026 
-24.965 ± 0.026 


52.88 ±0.16 
55.00 ±0.16 


10" 2 


KL' 
KL 


-11.082 ±0.054 
-2.302 ± 0.046 


80.309 ± 0.063 
79.838 ± 0.065 


-24.759 ± 0.015 
-24.759 ± 0.015 


44.465 ± 0.084 
52.566 ± 0.086 



Table 2: 



NLO Coefficient B(x) for cos\ = 0, y m i n = 10 4 


colour f actor 


region 


S(cos X = 0) 


r 2 


sing 


-410.18 ±0.16 




real 


408.36 ± 0.020 




fin 


-0.3744 ± 0.0019 




sum 


-2.18 ±0.26 


C F N C 


sing 


-68.818 ±0.071 




real(13) 


-120.62 ±0.12 




real (34) 


263.80 ±0.12 




fin(13) 


-0.3023 ± 0.0007 




fin(34) 


6.706 ± 0.016 




sum 


80.81 ±0.19 


C F N f 


sing 


-44.349 ± 0.018 




real 


19.227 ±0.037 




sum 


-25.116 ±0.041 


C 2 -C F N c /2 


(qq)int 


-0.0851 ± 0.0004 


sum 


sum 


53.52 ±0.32 
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Table 3: 



NLO Coefficient B(x) for cosx = 


Drain 


method 


cl 


C F N C 


C F N f 


sum 




KL' 
KL 


-2.67 ±0.68 
-2.58 ±0.68 


77.02 ± 0.60 
76.90 ± 0.60 


-25.388 ± 0.071 
-25.388 ± 0.071 


48.98 ±0.91 
48.95 ± 0.91 


10" 5 


KL' 
KL 


-2.03 ± 0.44 
-2.04 ± 0.44 


76.26 ± 0.35 
76.12 ±0.35 


-25.159 ±0.055 
-25.159 ±0.055 


49.06 ± 0.56 
48.92 ± 0.56 


10~ 4 


KL' 
KL 


-1.80 ±0.26 
-1.89 ±0.26 


74.39 ±0.19 
74.48 ±0.19 


-25.116 ±0.041 
-25.116 ±0.041 


47.47 ± 0.32 
47.49 ± 0.32 


10" 3 


KL' 
KL 


-1.73 ±0.13 
-1.86 ±0.13 


69.714 ± 0.086 
69.791 ± 0.087 


-24.965 ± 0.026 
-24.965 ± 0.026 


43.02 ±0.15 
42.97 ±0.15 


io- 2 


KL' 
KL 


-2.610 ±0.046 
-3.336 ± 0.049 


63.751 ± 0.041 
63.887 ± 0.046 


-24.759 ±0.015 
-24.759 ±0.015 


36.380 ± 0.063 
35.791 ± 0.068 



Table 4: 



NLO Coefficient B(x) fo r cos X = 


calculation 


B(cos x = 0) 


knw 


51.25 ±2.67 




52.39 ± 0.83 


GSW 


51.15 ±0.68 


GSW 


52.29 ±2.08 


KS& 


51.96 ±0.58 


KS& 


51.87 ±0.58 


KSW 


49.06 ± 0.56 


KS& 


48.92 ± 0.56 
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